Decays of rj and rj' within a chiral unitary approach 
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Various decays of and rj' are investigated within tiie framework of U(3) cliiral effective field 
theory in combination with a relativistic coupled-channels approach. Final state interactions 
are included by deriving s- and p-wave interaction kernels for meson-meson scattering from 
the chiral effective Lagrangian and iterating them in a Bethe-Salpeter equation. Very good 
agreement with experimental data is achieved. 
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1 Introduction 

Decays of 77 and rj' offer a possibility to study symmetries and symmetry breaking patterns in 
strong interactions. The isospin violating decays 77, rj' Svr, e.g., can only occur due to an 
isospin breaking quark mass difference m„ — or electromagnetic effects. While for most 
processes isospin violation of the strong interactions is masked by electromagnetic effects, these 
corrections are expected to be small for the three pion decays of rj and rj' (Sutherland's theorem) 
[ 1 ] . Neglecting electromagnetic corrections the decay amplitude is directly proportional to to„ — 
iTid- The decays 77,77' 77 and 77,77' 7r"'"7r~7, on the other hand, are phenomenological 
manifestations of the chiral anomaly of QCD and can provide important information on the 
chiral symmetry of the strong interactions. 

Moreover, the 77-77' system offers a testing ground for the role of both spontaneous and explicit 
chiral symmetry breaking, the latter one induced by the light quark masses. In the absence of 
77-77' mixing, 77 would be the pure member 773 of the octet of Goldstone bosons which arise due 
to spontaneous breakdown of chiral symmetry. 

Reactions involving the 77' might also provide insight into gluonic effects through the axial 
U(l) anomaly of QCD. The divergence of the singlet axial-vector current acquires an additional 
term with the gluonic field strength tensor that remains in the chiral limit of vanishing light quark 
masses. This term prevents the pseudoscalar singlet 770 from being a Goldstone boson which is 
phenomenologically manifested in its relatively large mass, 771^/ = 958 MeV. 

An appropriate theoretical framework to investigate low-energy hadronic physics is provided 
by chiral perturbation theory (ChPT) [2], the effective field theory of QCD. In ChPT Green's 
functions are expanded perturbatively in powers of Goldstone boson masses and small three- 
momenta. However, final state interactions, e.g., in 77 — > Stt have been shown to be substantial 
both in a complete one-loop calculation in SU(3) ChPT [3] and using extended Khuri-Treiman 
equations [4] . In 77' decays final state interactions are expected to be even more important due to 
larger phase space and the presence of nearby resonances. It is claimed, e.g., that the exchange 
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of the scalar resonance ao(980) dominates the decays 77' t^tttt [5] which has been confirmed 
both in a full one-loop calculation utilizing infrared regularization [6] and in a chiral unitary 
approach [7]. 

Also for the decays 77,77' and 77' — > 7r+7r~7 contributions from vector meson 

exchange will dominate the amplitude and unitarity effects should be implemented via final state 
interactions. The Dalitz decays 77, rj' 7*7 probe the charge distribution of 77 and 77', while the 
double Dalitz decays 77, 77' 7*7* may indicate whether double vector dominance is realized 
in nature. In particular the 77' decays are clearly beyond the perturbative framework of chiral 
perturbation theory (ChPT) and require utilization of non-perturbative tools which match onto 
the results from ChPT. 

In a series of papers [7-10] we have developed a theoretical framework which describes these 
77, 77' decays in a uniform manner. The approach is based on U(3) chiral effective field theory in 
combination with a relativistic coupled-channels approach. Final state interactions are included 
by deriving s- and p-wave interaction kernels for meson-meson scattering from the chiral effec- 
tive Lagrangian and iterating them in a Bethe-Salpeter equation (BSE). The illustration of the 
underlying idea of this method as well as a summary of the results for the hadronic decays is the 
subject of this presentation. 

2 Sketch of the approach 

Let us start with the hadronic decays 77, 77' Stt and 7/' — > rjinr [7, 10]. The underlying idea 
of the approach is that the initial particle, i.e. the rj or i]', decays into three mesons and that 
two out of these rescatter (elastically or inelastically) an arbitrary number of times, see Fig. ^ 
for illustration. All occurring vertices are derived from the effective Lagrangian and are thus 
constrained by chiral symmetry. Interactions of the third meson with the pair of rescattering 
mesons are neglected which turns out to be a good approximation, particularly for the decays 
77 — > Stt and 77' — > ijinr. 
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Fig. 1. Shown is a possible contribution to final state interactions in the decay rj —> n tv n . 

Such an infinite meson-meson rescattering process can be generated by application of the 
Bethe-Salpeter equation. To this aim, the partial wave interaction kernels for meson meson 
scattering, Ai, are derived from the chiral effective Lagrangian and iterated in a Bethe-Salpeter 
equation which generates the propagator for two interacting particles in a covariant fashion. For 
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each partial wave I the matrix-valued solution T; of the BSE with coupled channels and on-shell 
interaction kernels is given by 

T,-' = Ai^ + G , (1) 

where the diagonal matrix G collects the scalar loop integrals of the different two-meson chan- 
nels. The solution of the BSE satisfies exact unitarity for two-particle scattering and gener- 
ates resonances dynamically by an infinite string of meson-meson rescattering processes, see 
also [11]. In the investigations discussed here we have restricted ourselves to s- and p-wave 
interaction kernels. 

This approach can be extended to the anomalous decays rj, 77' — » [8] and 77, 77' 

7r+7r^7 [9] in a straightforward manner In ?/, 77' — > 7(*)'y(*), e.g., the incoming pseudoscalar 
meson P can directly decay via a vertex of either the Wess-Zumino-Witten Lagrangian, which 
describes the chiral anomaly, or the unnatural-parity Lagrangian into one of the following three 
channels: two photons, a photon and two pseudoscalar mesons, or four pseudoscalar mesons. 
Pairs of mesons can then rescatter an arbitrary number of times before they eventually couple to 
a photon, see Fig.|2]for illustration. In order to implement the non-perturbative summation of 
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Fig. 2. Sample rescattering process in the decay P -y(*)'y(*). 

loop graphs covered by the BSE in the decay processes under consideration, we treat the BSE 
T-matrix as an effective vertex for meson-meson scattering. The anomalous decays into 7r+7r~7 
are treated in the same way [9]. 

3 Results 

We now turn to the discussion of the numerical results, where for brevity we restrict ourselves to 
the hadronic decays [10]. Results for the anomalous decays can be found in [8,9]. The results 
of our calculation for the hadronic decays are obtained from a combined analysis of the decay 
widths, branching ratios, and slope parameters of the considered decays as well as phase shifts in 
meson-meson scattering. The widths of 77 ^ Stt and 77' — > iitttt have been measured roughly at 
the 10 % precision level, while for r/' Svr*' the experimental uncertainty is considerably larger 
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and only an upper limit exists for r(77' tt+tt 7r°) [12]. Moreover, some of these decay widths 
are constrained by the well-measured branching ratios 

'^^ r(7/ ^ tt+tt-ttO) ' ""^ r(7/ -> ?77r"7rn) ' ^ ' 

The squared absolute value of the decay amplitude, is expanded for rj' rj-wK 

and the charged decay modes of -q, rj' 3tt in terms of the usual Dalitz variables x, y, as 

\Aix,y)f = \N\^[l + ay + by^ + cx^ + dy^ + ---] , (3) 

while for the decays into three identical particles Bose symmetry dictates the form 

\A{x,y)\'^\N'\'[l + g(y' + x') + ■■■]. (4) 

For the Dalitz plot parameters a, b, c of ?; ^ tt+tt^tt'^ the experimental situation is not without 
controversy. We employ the numbers of [13], since it is the most recent published measurement 
and the results appear to be consistent with the bulk of the other experiments listed by the Par- 
ticle Data Group [12]. They differ somewhat from the new preliminary results of the KLOE 
Collaboration [14] that has found a non-zero value for the third order parameter d which was not 
included in previous experimental parametrizations. Very recently the Dalitz plot parameters of 
rj' r^TT+TT^ have been determined with high statistics by the VES experiment [15]. In this 
presentation, we will discuss the fits obtained from the experimental Dalitz parameters provided 
by the Particle Data Group [12]. Note that the slope parameters of 77' Sir have not yet been 
determined experimentally, but such a measurement is intended at WASA@COSY [16]. 

From the unitarized partial-wave T-matrix one may also derive the phase shifts in meson- 
meson scattering. Hence, our approach is further constrained by the experimental phase shifts 
for TTTT TTTT, KK Scattering. 

The coupled-channels framework entails several chiral couplings of the effective Lagrangian 
up to fourth order which must be fit to experimental data. By fitting to all available (published) 
data sets of the investigated hadronic rj, ?/ decays and phase shifts an overall function is 
calculated. We observe four different classes of fits which are all in very good agreement with 
the currently available (published) data on hadronic decays, but differ in the description of the 
decays rj' 37r where experimental constraints are scarce. The errors which we specify in the 
following for all parameters and observables reflect the deviations which arise when we allow 
for values which are at most 15 % larger than the minimum value. Although this choice is 
somewhat arbitrary, it illustrates how variation of the function in parameter space affects the 
results. The results for 77 Sir are shown in Tabled 

Only sparse experimental information exists on the decays of rj' into three pions. The exper- 
imental decay width of rj' — > 37r° is [12] 

which is nicely met within our approach: 

p(theo)(^/ ^ 3^0) ^ (330 ± 33)ey (6) 
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r„^Q.,o (eV) 


r„ ^+^-^0 (eV) 


' ± 




theo. 


422 ±13 


290 ±8 


1.46 ±0.02 


exp. 


419 ±27 


292 ±21 


1.44 ±0.04 




a 


b 


c 


9 


theo. 


-1.20±0.07 


0.28 ±0.05 


0.05 ±0.02 


-0.062 ±0.006 


exp. 


-1.22±0.07 


0.22±0.11 


0.06 ±0.02 


-0.062 ±0.008 



Tab. 1. Results for the partial decay widths of — > Sir, the branching ratio ri, and the Dalitz plot parameters 
compared to experimental data from [12] and [13]. 



Forthedecayint07r+7r~7r° onlyaweak experimental upper limitexists, r^'^''P^(77' — > 7r+7r~7r'^) < 
lOkeV [12]. As already mentioned above, we observe four different clusters of fits which 
lead to equally good agreement with the available data but yield different predictions for the 
7]' 7r+7r~7r° decay width and the urmieasured rj' — » Stt Dalitz parameters, see [10] for de- 
tails. The Dalitz plot distributions of these decays and the decay width r]' 7r+7r~7r° pose 
therefore tight constraints for our approach and must be compared with future experiments at the 
WASA@COSY facihty. 

For Tn'^n2n we obtain the results 

^Oheo) ^ f^j^ ^js^^ ^(exp) ^ ^ x 10"^ [12]. (7) 





77'— >777r+7r~ 


a 


b 


c 


theo. 


81 ±4 keV 


-0.116 ±0.024 


0.000 ±0.019 


0.016 ±0.035 


exp. 


89±llkeV 


-0.16 ±0.06 







Tab. 2. Results for the partial decay width of rj' rjir'^'K and the DaUtz plot parameters compared to 
experimental data from [12]. 







a 


b 


c 


theo. 


46±3keV 


-0.122 ±0.025 


0.003 ±0.018 


0.019 ±0.039 


exp. 


42 ± 6 keV 


-0.116±0.026 


0.003 ±0.017 


0.00 ±0.03 



Tab. 3. Results for the partial decay width of r;' rin°7r° and the Dalitz plot parameters compared to 
experimental data from [12]. 

In the isospin limit, m„ = ma, the decay width r{r]' 7777+77") would be exactly given 
by 2 r{ri' r^Tr^Tr^), due to the symmetry factor for identical particles in the latter process. If, 
however, one is interested in the isospin-breaking contributions in the amplitude of 77' rjmr, 
one ought to disentangle it from phase space effects which are caused by the different masses of 
charged and neutral pions. With an isospin-symmetric decay amplitude, but physical masses in 
the phase space factors, we find a ratio 
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which is smaller than 2 and compares to rs = 1.77 ± 0.02 when isospin-breaking is taken into 
account in the amplitude. (For comparison, if the amplitude is set constant and the physical pion 
masses are employed in the phase space integrals, the ratio is given by rs = 1.77.) We may thus 
conclude that within our approach isospin-breaking corrections in the rj' rjinr decay amplitude 
are tiny. The branching ratio ra has not been measured directly. If, however, we calculate the ratio 
of fractions Fj/rtotai for these two decay modes using the numbers and correlation coefficients 
pubhshed by the Particle Data Group [12], we arrive at 

rf'"'^^ = 2.12 ± 0.19 (9) 

by means of standard error propagation. Such a large branching ratio would indicate significant 
isospin- violating contributions in the amplitude. But the experimental uncertainties are sizable 
and should be reduced by the upcoming experiments with WAS A at COSY [16] and at MAMI- 
C[17]. 

To summarize the results for the hadronic decays of ij and rj' in our approach, we obtain very 
good agreement with currently available data on the decay widths and spectral shapes. In fact, 
we observe four different classes of fits which describe these data equally well, but differ in their 
predictions for yet unmeasured quantities such as the rj' 7r+7r^7r° decay width (for which 
there exists only a weak upper limit) and the Dalitz slope parameters of rj' — > Stt. The results 
obtained may be tested in future experiments foreseen at WASA@COSY and MAMI-C. The 
hadronic decays considered here along with phase shifts in meson-meson scattering pose there- 
fore tight constraints on the approach and will allow to determine the couplings of the effective 
Lagrangian up to fourth chiral order It is important to stress that the values of the parameters 
obtained from the fit are in general not the same as in the framework of ChPT which can be 
traced back to the absorption of loops into the coefficients and higher order effects not included 
in the perturbative framework. 

An intriguing feature of the fits is that they accommodate the large negative slope parameter 
g of the decay r] — > 37r° measured by the Crystal Ball Collaboration [18] which could not be met 
by previous theoretical investigations [4,7]. This value must, however, be confronted with the 
more recent but yet prehminary g value of the KLOE Collaboration [14]. If we replace the PDG 
data by the KLOE Dalitz parameters of both the charged and neutral — > Stt decay, we do not 
achieve a good overall fit. It appears that the slope parameters of both ry ^ Stt decays cannot 
be fitted simultaneously. In addition, fitting to the KLOE data destroys the agreement with the 
experimental branching ratio of both decays which is known to high precision. To this end, we 
have illustrated that utilizing the AI = 1 selection rule which relates both decays and taking the 
KLOE parametrization of the charged decay as input leads in a model-independent way to a g 
value not consistent with the KLOE g result. 

The importance of the various two-particle channels with different isospin and angular mo- 
mentum has been examined as well. For the rj ^ Sn decays we find that the major contribution 
is given by tttt rescattering in the s-wave / = channel, while the J = 1, 2 channels interfere 
destructively with the former. The p-wave contribution in the charged decay is tiny, since avail- 
able phase space is small and the C-odd channels related to the p(770) resonance do not occur. 
For T]' — > Tr~^TT~TT^, on the other hand, phase space is considerably larger, and the size of the 
p-wave contributions ranges from 10 % to 50 % depending on the cluster of fits. 

For the decays rj' tjittt we find that the s-wave 1 = 1 channels dominate for two classes of 
fits which would confirm the importance of the nearby ao (980) resonance as claimed by previous 
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investigations [5-7]. But the other two clusters are dominated by the / = channels. These two 
scenarios can be distinguished by their predictions for the ij' Stt decays. Thus, a precise 
measurement of 77' Stt decay parameters can also help to clarify the importance of ao(980) or 
/o(980) resonance contributions to the dominant decay mode of the 77' into rjinr. 

Acknowledgement: This work was supported in part by DFG, SFB/TR-16 "Subnuclear Struc- 
ture of Matter", and Forschungszentrum Jiilich. 
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